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1.  Introduction 

Since  plasma-boundary  physics  encompasses  some  of  the  most 
important  unresolved  issues  in  future  energy  production  for  fusion 
reactors,  there  is  a  strong  interest  in  the  fusion  community  for  bet¬ 
ter  understanding  and  characterization  of  plasma-surface  interac¬ 
tions  (PSI).  Chemical  and  physical  sputtering  cause  the  erosion  of 
the  limiter/divertor  plates  and  vacuum-vessel  walls,  whether  these 
are  made  of  C,  Be,  W,  or  some  combination  (i.e.  mixed  materials), 
and  degrade  fusion  performance  by  diluting  the  fusion  fuel  and 
excessively  cooling  the  core.  Hydrocarbon  re-deposition  onto  plas¬ 
ma-facing  components  can  lead  to  long  term  accumulation  of  large 
in-vessel  tritium  inventories. 

A  major  challenge  in  the  production  of  D-T  fusion  power  is  the 
development  of  materials  for  the  first  wall  and  internal  compo¬ 
nents.  The  choice  of  wall  material  has  profound  effects  on  confine¬ 
ment  of  fusion-grade  plasmas.  Although  carbon-based  materials 
have  superior  thermo-mechanical  properties,  they  could  trap  high 
levels  of  tritium  by  co-deposition  with  eroded  carbon  and  thereby 
severely  constrain  plasma  operations.  Thus,  a  mix  of  several  differ¬ 
ent  plasma-facing  materials  is  now  proposed  in  ITER  to  optimize 
the  requirements  of  areas  with  different  power  and  particle  flux 
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characteristics.  The  slow  rate  of  progress  in  the  area  of  tritium  re¬ 
moval,  together  with  favorable  results  from  divertor  tokamaks 
with  high  atomic  number  (e.g.  tungsten)  walls,  suggest  that  inter¬ 
est  in  all-metal  machines  will  increase. 

The  key  issues  for  surface  experiments  on  carbon-based  materi¬ 
als  are  erosion,  reflection,  impurity  transport  in  the  plasma,  depo¬ 
sition,  T  uptake  and  removal.  There  is  a  need  to  determine  the 
composition  of  eroded  species  such  as  hydrocarbon  molecules 
and  radicals,  their  rovibrational  state  and  energy  spectra,  and  their 
sticking  coefficients  to  surfaces  as  a  function  of  energy.  An  expan¬ 
sion  of  the  available  erosion  database  towards  very  low  energies 
(5  eV)  is  needed  in  order  to  narrow  the  gap  in  erosion  data  between 
energetic  hydrogen  ions  and  thermal  atomic  hydrogen,  both  for 
pure  carbon  as  well  as  mixed  materials  systems.  Properly  bench- 
marked  molecular  dynamics  (MD)  simulations  could  provide  com¬ 
prehensive  databases  for  boundary  plasma  modeling,  providing 
details  often  not  accessible  by  experiment. 

Chemical  sputtering  is  a  process  where  bombardment  by  atoms 
or  molecules  induces  a  chemical  reaction  which  produces  a  particle 
that  is  weakly  bound  to  the  surface  and  hence  is  easily  desorbed 
into  the  gas  phase  [1].  For  carbon  surfaces,  it  is  hypothesized  that 
incident  particles  break  bonds  within  the  collision  cascade  [2].  The 
broken  bonds  are  rapidly  passivated  by  the  abundant  flux  of  atom¬ 
ic  hydrogen  from  the  hydrogen  fusion  plasma  environment. 
This  leads  to  the  formation  of  stable  hydrocarbon  molecules 
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underneath  the  surface,  which  diffuse  to  the  surface  and  desorb 
thermally. 

In  our  MD  simulations  the  classical  motion  of  a  collection  of 
atoms  is  followed  using  forces  obtained  from  predefined  analytic 
ground  electronic  state  potentials,  which  for  hydrocarbon  systems 
have  reached  maturity  (empirical  bond-order  potentials  of  the 
Tersoff-Brenner  type,  such  as  REBO  [3,4]  and  AIREBO  [5])  and  their 
predictive  strength  is  improving.  The  parameters  of  the  phenome¬ 
nological  potential  functions  are  taken  either  from  parameter-free 
quantum-mechanical  calculations  or  experimental  data.  Dynamic 
electronic  effects  are  not  included  in  classical  MD.  Non-local  ef¬ 
fects,  which  take  into  consideration  the  bonding  structure  further 
away  than  the  nearest  neighboring  atoms,  may  be  included  in 
the  potential  to  improve  the  accuracy  of  the  calculations.  This  is 
especially  important  for  carbon,  where  heterogeneous  bonding 
configurations  arise  from  the  presence  of  hybridized  bonding. 
The  REBO  potential  has  a  reasonably  realistic  description  of  pure 
carbon  and  hydrocarbon  molecular  structures,  as  well  as  dynamic 
effects,  such  as  bond  forming  and  breaking.  Hence,  it  has  been 
adopted  by  most  MD  studies  of  sputtering  [6-9]  including  our 
own  [10-15].  The  AIREBO  potential  overcomes  some  limitations 
of  the  REBO  potential,  but  requires  significantly  faster  computer 
resources,  due  to  the  longer  range  of  the  potential. 


CxDy  / 


Fig.  1.  (a)  Chemical  sputtering  upon  impact  of  D  and  D2  projectiles  on  deuterated 
amorphous  carbon  and  (b)  filamentous  low-density  structures  at  the  interface  after 
prolonged  bombardment  with  20  eV  D. 


The  numerical  cell  must  be  larger  than  the  characteristic  region 
of  atom  impact  induced  perturbation,  to  secure  scattering-free  cell 
boundaries.  In  the  range  we  have  studied  [10-15],  below  30  eV  im¬ 
pact  energy,  a  realistic  cell  of  deuterated  carbon  consists  of  a  cubic 
random  C:D  network  of  several  thousand  atoms,  with  edge  dimen¬ 
sion  of  2.5  nm  (Eig.  1(a))  and  two-dimensional  periodic  boundary 
conditions  in  the  planes  perpendicular  to  the  surface  plane.  Satu¬ 
rated  surfaces  are  prepared  by  cumulative  bombardment  of  several 
thousand  projectiles  during  creation  times  up  to  10  ns  until  a  sat¬ 
uration  regime  is  reached  [10,15].  Cumulative  bombardment  not 
only  erodes  the  surface,  but  also  creates  surface  swelling  producing 
a  lower-density  interface  that  is  inhomogeneous  and  contains  fila¬ 
mentous  structures  [14]  of  hydrocarbon  chains  (Eig.  1(b)),  which 
play  an  important  role  in  the  PSl  dynamics.  Bombardment  also 
changes  the  D/C  ratio  close  to  the  surface  interface  (extending  to 
increasing  depths  with  increasing  impact  energy),  where  it  reaches 
values  ^1,  as  seen  in  Eig.  2.  Sputtering  yields  typically  reach  a 
nearly  stationary  regime  as  a  function  of  fluence.  However,  this 
stationarity  can  never  be  a  true  steady-state  since  the  finite  simu¬ 
lation  cell  is  continually  changing  due  to  erosion.  In  any  case, 
chemical  sputtering  is  strongly  dependent  on  the  accumulation 
of  hydrogen  in  the  surface  and,  consequently,  depends  on  the  con¬ 
centration  of  carbon  atoms  with  sp/sp^/sp^  hybridization  states, 
which  also  changes  with  fluence.  This  explains  the  change  in  the 
total  carbon  yield  as  a  function  of  the  deuterium  fluence,  shown 
in  Eig.  3,  which  is  also  seen  experimentally  [16].  Most  of  our  sim¬ 
ulations  were  done  with  surfaces  prepared  in  the  range  of  1000- 
2000  impacts  (at  or  beyond  the  point  where  the  sputter  yield  has 
reached  an  approximate  stead  state  in  Eig.  3). 

After  creating  the  surfaces  by  cumulative  bombardment,  sput¬ 
tering  yields  are  subsequently  calculated  by  independent  random 
impacts  of  several  thousand  projectiles  on  the  saturated  surfaces 
[10].  At  the  lowest  energies  the  time  duration  of  the  collisional  cas¬ 
cade  for  each  impact  is  less  than  5  ps,  but  approaches  30  ps  for 
30  eV  deuterium  impacts  [15].  Calculations  are  performed  using 
2000-4000  independent  trajectories  for  5-30  ps.  Such  level  of  sta¬ 
tistics  is  essential  for  meaningful  calculation  of  small  yields.  Our 
recent  calculations  [10-15]  of  impact  of  deuterium  atoms  and  mol¬ 
ecules  have  shown  good  agreement  with  the  experimental  results 
[17]  for  sputtering  of  methane  and  acetylene.  This  good  agreement 
for  selected  hydrocarbon  yields  is  the  result  of  our  computational 
methodology:  we  mimic  the  cumulative  saturation  conditions  of 


Fig.  2.  Depth  profile  for  the  density  of  carbon  (open  symbols)  and  deuterium  (solid 
symbols)  in  the  simulation  cell,  both  for  the  initial  surface  (black  squares)  and  for 
the  surface  after  600  impacts  (red  circles).  Initially,  the  cell  was  created  from  a 
random  distribution  with  a  D/C  ratio  of  0.4,  homogenized  by  a  heating  and 
annealing.  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the 
reader  is  referred  to  the  web  version  of  this  article.) 
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Fig.  3.  Evolution  of  carbon  yield  with  increasing  D  fluence  in  the  process  of  surface 
preparation. 


the  beam-surface  experiment  as  closely  as  possible,  preparing  the 
surface  self-consistently  by  bombardment  with  appropriate  parti¬ 
cles  at  appropriate  energies,  until  a  quasi  steady-state  regime  in 
the  sputtering  yields  is  reached  [10].  This  results  in  hundreds  of 
prepared  surfaces,  for  various  energies  and  impact  particles,  which 
are  then  used  in  the  computational  simulations  of  chemical  sput¬ 
tering,  repeated  for  each  case  thousands  of  times,  and  averaged 
over  a  number  of  surfaces,  obtained  at  various  fluences.  We  have 
calculated  particle,  energy  and  angular  spectra  of  the  sputtered 
hydrocarbons,  reflection  coefficients,  penetration  coefficients,  dis¬ 
sociation  probabilities  of  impacting  D2,  and  D2  production  at  the 
surface  with  impact  of  D  [15].  We  have  also  shown  that  the  sput¬ 
tering  yields  of  hydrocarbons  strongly  depend  on  the  initial  vibra¬ 
tional  states  of  impacting  D2  molecules  [11]. 

In  addition  to  particle  beams,  well-controlled  plasma  sources 
can  also  be  used  for  experimental  characterization  of  sputtering 
processes.  Experiments  of  this  kind  were  performed  recently  [18] 
for  the  accommodation  coefficients  of  H2  interacting  with  poly¬ 
crystalline  carbon  surfaces.  We  have  therefore  calculated  the  sput¬ 
tering  yields  and  reflection  coefficients  using  a  distribution  of 
impacting  particles  that  mimics  a  plasma  environment,  which 
can  be  directly  validated  by  experiment. 

In  the  following  we  present  a  brief  review  of  typical  results  ob¬ 
tained  in  our  simulations  of  surface  irradiation  with  particle  beams 
and  new  results  with  the  plasma-irradiated  surfaces. 

2.  Results 

Most  of  our  results  were  obtained  under  the  assumption  of  a 
monoenergetic  beam  irradiation,  incident  perpendicular  to  the  sur¬ 
face,  which  is  kept  at  temperature  of  300  K  using  a  Langevin  ther¬ 
mostat  (for  details  see  [10-15,17]).  We  have  studied  the  least  well 
known,  intermediate-to-low  range  of  impact  energies  (1-30  eV/D), 
with  D  and  D2  projectiles  impacting  deuterated  amorphous  carbon 
surfaces. 

2.1.  Chemical  sputtering  and  surface  microstructure 

Assuming  the  bond  breaking  as  the  first  step  in  the  chemical 
sputtering  process,  hydrocarbon  moieties  R-CDx  at  the  interface 


Impact  energy  (eV/D) 

Fig.  4.  (a)  Counts  of  R-CDx  moieties  in  the  simulation  cell  as  a  function  of  impact 
energy,  integrated  over  the  portion  of  the  simulation  cell  with  z  >  0,  for  a  range  of 
impact  energies.  The  counts  represent  averages  across  six  surfaces  [10].  (b)  Sputter 
yield  of  CD^  hydrocarbons. 


which  are  created  at  the  terminal-ends  of  the  branched  network 
structure  are  expected  to  play  a  dominant  role  in  the  sputtering 
of  the  simplest  hydrocarbons,  like  CD,  CD2  and  CD3.  Furthermore, 
the  sputtering  yield  dependence  of  these  hydrocarbons  on  impact 
energy  should  be  proportional  to  the  corresponding  densities  of  sp, 
sp^  and  sp^  terminal  carbons,  respectively,  generated  in  the  surface 
at  the  various  sputtering  yields.  Our  simulations  show  striking 
similarities  between  the  shapes  of  the  relevant  yields  and  the  sp^ 
densities  as  functions  of  the  impact  energy/D  (Fig.  4),  confirming 
that  the  leading  mechanism  for  the  CDx  sputtering  is  the  bond 
breaking  of  terminal-end  R-CD^  moieties  [10]. 

2.2.  Vibrational  excitation  of  the  impinging  molecules  [11] 

The  impacting  particle  is  most  efficient  in  chemical  sputtering 
at  the  end  of  its  thermalization  cascade.  Penetration  of  an  impact¬ 
ing  molecule  (like  D2)  into  a  hydrogenated  carbon  surface  might  be 
suppressed  due  to  the  molecule  physical  dimensions.  Thus,  the  dis¬ 
sociation  of  the  molecule  is  a  predecessor  for  the  successful  chem¬ 
ical  sputtering.  However,  depending  on  the  vibrational  excitation 
of  the  D2  molecule,  a  fraction  of  the  impact  energy  must  be  ex¬ 
pended  for  the  dissociation,  which  depends  on  how  highly  the 
molecule  is  excited.  Therefore,  the  sputtering  yield  increases  with 
increasing  vibrational  excitation,  as  shown  in  Fig.  5,  as  function  of 
the  classical  vibrational  amplitude  of  the  D2  molecule.  The  maxi¬ 
mum  is  reached  at  the  end  of  the  correlation  distance  of  the  short 
range  REBO  potential,  about  2  A.  Beyond  the  maximum,  when  the 
atoms  are  not  directly  correlated,  the  sputtering  yield  decreases 
slightly  since  the  interaction  between  the  two  impact  cascades  de¬ 
creases  with  distance.  When  the  distance  between  the  impacting 
atoms  becomes  large  enough,  the  sputtering  yield  tends  to  twice 
the  single  D-atom  yield,  also  shown  in  Fig.  5. 

2.3.  Validation  of  the  MD  simulation:  comparison  with  the  beam 
experiments 

By  mimicking  experiments  through  preparation  of  the  surfaces 
self-consistently  in  energy  and  particle  type  by  cumulative  bom¬ 
bardment,  very  good  agreement  was  obtained  between  the  MD 
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Vibrational  amplitude  (Angstroms) 


Fig.  5.  Calculated  sputtering  yields  per  deuterium  atom  of  acetylene  (C2D2,  open 
triangles)  and  total  carbon  (filled  squares)  for  a  D2  impact  energy  of  10  eV/D,  as  a 
function  of  the  initial  internal  vibrational  state  of  D2,  expressed  through  the 
maximum  “bond”  length  (classical  vibration  amplitude),  1  [11].  The  yields  for 
impact  of  non-interacting  D  are  multiplied  by  a  factor  of  two  and  are  shown  by  the 
rightmost  horizontal  lines. 

simulation  results  for  the  methyl  +  methane  sputtering  yield  and 
measurements  of  the  ejected  methane  yield  [10,17]  (Fig.  6).  We 
note  that  the  experimental  results  in  Fig.  6  correspond  to  incident 
Dj,  which  could  be  mimicked  by  two  D  atoms  (see  Fig.  5)  if  the  D2, 
obtained  by  neutralization  of  the  Dj  above  the  surface,  is  in  a 
highly  vibrationally  excited  state.  It  is  not  currently  clear  what  is 
the  role  of  the  post-sputtering  processes  in  the  experimental 
chamber,  nor  whether  our  simulation  contains  a  complete  descrip¬ 
tion  of  methane  production.  However,  it  is  of  no  concern  that  the 
measurements  were  done  with  Dj  or  D\  and  simulations  with 
D2  (or  with  D):  The  neutralization  of  Dj  is  expected  to  happen 
through  resonant  charge  transfer  to  vibrationally  excited  D2  [11], 
and  neutralization  of  goes  through  resonant  charge  transfer  to 
D,  which  are  very  fast  processes  in  both  cases. 

2.4.  Mass  spectra  of  the  ejected  hydrocarbons  [10,12,15] 

The  preparation  of  the  surfaces  by  cumulative  bombardment  at 
each  energy  and  particle  type  turned  out  to  be  essential  for  sput¬ 
tering  of  complex  hydrocarbons  [10].  A  “virgin”  surface  has  the 
tendency  to  eject  only  carbon  atoms  and  diatoms,  and  the  very 
simplest  hydrocarbons  (like  CH).  Fig.  7  shows  the  mass  spectra  of 
hydrocarbons,  ejected  in  MD  simulations  by  vibrationally  excited 


Fig.  6.  Comparison  of  the  MD  result  for  sputtering  of  CD3  +  CD4  by  D  impact  with 
the  results  of  beam  experiments  for  D2  impact  [19]. 


Mg  (amu) 


Fig.  7.  Yields  of  hydrocarbons,  sputtered  by  D2  impact  on  a-C:D,  as  functions  of 
mass  of  the  ejected  molecules,  for  various  impact  energies  [10,12]. 

D2.  The  groups  QDy  are  well  separated  in  x.  With  increasing  impact 
energy,  more  and  more  complex  hydrocarbons  arise.  While  at  the 
lowest  energies  mostly  x  =  1  hydrocarbons  are  sputtered,  at 
30  eV/D  hydrocarbons  with  up  to  x  =  4  are  visible  in  the  mass  spec¬ 
trum.  The  distribution  of  the  ejected  hydrocarbons  is  not  uniform 
withy  for  a  given  x  but  peaks  at  progressively  less  saturated  hydro¬ 
carbons  with  increasing  mass.  For  x  =  1  the  biggest  yield  is  of 
methyl,  CD3,  for  x  =  2  it  is  for  acetylene,  C2D2,  for  x  =  3  the  maxi¬ 
mum  is  at  C3D4,  and  so  on. 

2.5.  Energy  spectra  of  the  hydrocarbons  [10,12,15] 

Energies  of  the  ejected  hydrocarbons  as  function  of  their  mass, 
for  various  values  of  the  impact  D2  energy,  are  shown  in  Fig.  8. 
Although  the  dispersion  of  energies  seems  to  decrease  with  mass, 
it  is  rather  counterintuitive  that  both  light  and  heavy  ejected  par¬ 
ticles  have  similar  average  energies,  irrespective  of  the  impact  en¬ 
ergy.  These  average  values  are  between  0.5  and  1  eV,  and  are  too 
energetic  to  arise  thermally  from  the  surface  at  300  K.  An  explana¬ 
tion  is  that  sputtering  is  a  kinetic  process,  i.e.  a  hydrocarbon  is 
ejected  by  a  chemical  processes  at  the  end  of  the  collisional  cas¬ 
cade  [10].  The  sputtered  particle  is  thermalized  with  the  locally 


Fig.  8.  Energy  spectra  of  ejected  hydrocarbons  as  functions  of  hydrocarbon  mass 
[10,12].  Each  point  indicates  an  individual  ejecta.  The  solid  line  represents  the 
average  at  each  mass  number. 
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Fig.  9.  Yields  of  D,  D  in  D2  and  D  in  hydrocarbons,  and  their  sum  (total  D)  for 
various  impact  energies  of  atomic  D  at  normal  incidence  [15]. 


elevated  temperature,  which  is  confirmed  by  the  Maxwell-Boltz- 
mann  distribution  of  the  ejected  particles  [10,15]. 

2.6.  Balance  of  various  channels  of  hydrogen  ejection  [10,15] 

Deuterium  is  ejected  in  the  form  of  reflected  (and  sputtered) 
atoms,  as  well  as  in  the  form  of  sputtered  D2  and  hydrocarbons 
(Fig.  9).  At  impact  energies  in  the  range  4-5  eV,  there  is  a  change 
in  trends  of  the  curves  in  Fig.  9  which  corresponds  to  the  covalent 
bond  dissociation  energy  in  a-C:D.  Above  5  eV,  the  reflection  coef¬ 
ficient  continues  to  decrease,  although  much  more  slowly.  A  larger 
fraction  of  the  D-atoms  is  retained  in  the  substrate  at  higher  ener¬ 
gies.  The  D2  yield  is  largest  at  5-10  eV  and  decreases  gradually  as 
the  impact  energy  increases  above  10  eV.  The  continued  decrease 
in  D  reflection  is  due  both  to  the  increased  D  retention  (implanta¬ 
tion)  as  well  as  growth  in  the  yield  of  hydrocarbons  that  carry  off 
multiple  D-atoms. 

In  presence  of  erosion,  one  might  even  expect  that  total  ejection 
yield  of  hydrogen  (in  all  channels)  could  be  even  larger  than  1.  This 
does  not  happen  here,  indicating  that  the  surface  is  accreting 
hydrogen  faster  than  it  is  eroding  at  these  fluences.  There  is  erosion 
of  carbon,  which  increases  with  impact  energy,  but  the  yield  of 
hydrocarbons  stays  low  (as  is  well  known). 

2.7.  Angular  distribution  of  ejecta  [10,15] 

Because  the  target  is  amorphous  and  we  only  consider  D  im¬ 
pacts  normal  to  the  surface,  the  angular  distributions  dNjdQp  are 


independent  of  the  azimuthal  angle  in  the  plane  of  the  surface 
interface.  Therefore,  for  the  center  of  mass  (c.m.)  momentum  P  of 
the  ejected  particle,  only  the  dependence  on  the  polar  angle  0p  with 
respect  to  the  surface  normal,  is  of  interest. 

The  angular  distributions  for  sputtered  D2  and  hydrocarbons, 
for  various  impact  energies  of  D  are  presented  in  Fig.  10  [15].  For 
physical  sputtering  and  normal  incidence  on  an  amorphous  solid 
the  angular  distributions  of  the  ejected  momenta  are  expected  to 
follow  dNjdQp  ~  cosOp.  The  present  distributions  for  chemical  sput¬ 
tering  approximately  mimic  this  behavior,  particularly  at  high 
energies  although  the  curves  are  somewhat  more  sharply  peaked 
around  ejection  normal  to  the  surface  than  cosOp. 

2.8.  Rovibrational  spectra  of  sputtered  molecules  [15] 

Rotational,  vibrational  and  rovibrational  energies  of  ejected  D2 
and  hydrocarbons  are  presented  in  Fig.  11  as  functions  of  the 
impact  energy  of  D.  The  features  of  the  curves  in  the  range  of 


Fig.  11.  Average  ejection  energy,  including  internal  energy  (vibrational  and 
rotational)  of  ejected  molecules  by  D  impact,  for  ejected:  (a)  D2  and  (b)  hydrocar¬ 
bons  [15]. 


a  D  impact 


b  D  impact 


Fig.  10.  Angular  distributions  of  the  c.m.  momentum  of  ejected  molecules  upon  D  impact  [15]  as  a  function  of  the  polar  angle  0p  with  respect  to  the  surface  normal.  The  thick 
dashed  lines  represent  an  angular  distribution  dNjdQp  ~  cosOp.  The  distributions  are  normalized  such  that  the  maximum  values  of  dNjdQ  lie  on  the  unit  circle. 
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Fig.  12.  Rovibrational  energy  distribution  of  emitted  hydrocarbons  following  20  eV 
D  impact.  The  line  is  a  Boltzmann  distribution  ^expi-EjEu)- 


5-10  eV  impact  energy  reflect  the  C-D  and  C-C  bond  breaking  en¬ 
ergy  in  4-5  eV  range.  The  translational  energies  of  ejected  mole¬ 
cules  are  in  range  of  0.5-1  eV.  In  the  same  interval  are 
rovibrational  energies  of  D2  (Fig.  11(a)),  approximately  equiparti- 
tioned  between  rotational  and  vibrational  energies,  especially  at 
the  higher  end  of  impact  energies.  In  turn,  rovibrational  energies 
of  hydrocarbons  are  significantly  higher  than  translational  ener¬ 
gies,  about  2  eV,  mostly  due  to  the  increased  number  of  vibrational 
modes  of  the  complex  molecules.  Of  course,  the  corresponding 
temperatures  per  degree  of  freedom  are  closer  to  the  rotational 
and  translational  temperatures,  in  the  range  of  5000  K.  As  can  be 
seen  from  Fig.  12,  the  motion  of  ejecta  is  at  least  approximately 
thermalized  at  an  elevated  temperature,  similarly  to  the  transla¬ 
tional  motion  of  the  ejected  molecules  [10]. 

2.9.  “Plasma”  irradiation 

Simulations  have  shown  that  it  is  important  to  properly  mimic 
the  experimental  conditions  and  history.  For  a  surface  exposed  to  a 
plasma  this  includes  bombarding  the  surface  with  random  angles 
and  energy  energies  of  impact.  We  mimic  a  plasma  by  randomly 
sampling  the  impinging  energies  of  hydrogen,  deuterium  and  tri¬ 
tium  atoms  from  a  Maxwell  distribution,  and  the  impact  angles 
from  a  uniform  distribution.  Here,  we  show  preliminary  results, 
which  acquired  enough  statistical  weight  only  for  atom  reflection 
and  sputtering  of  hydrogen  molecules.  Each  surface  used  as  a  tar¬ 
get  was  prepared  self-consistently  by  randomized  cumulative  irra¬ 
diation  with  a  large  number  (1000-2000)  of  like  atoms,  at  a  desired 
plasma  temperature.  Results  for  ejection  of  impact-like  particles 
and  sputtering  of  impact-like  diatomic  molecules  are  shown  in 
Fig.  13,  for  a  range  of  sampled  temperatures  10,000-200,000  K 
(0.86-17  eV).  It  is  interesting  to  note  that  the  reflection  yield  stays 
(independent  of  the  isotopic  composition)  at  the  highest  tempera¬ 
tures  in  the  range  0.5-0.6  after  the  initial  decline,  somewhat  higher 
than  that  for  the  normal-beam  irradiation  (Fig.  9)  at  the  highest 
energies.  This  is  due  to  the  averaging  over  impinging  angles  and 
energies  that  takes  place  for  a  plasma-surface  interaction.  For  the 
same  reason  sputtering  yields  of  hydrogen  molecules  stay  below 
0.1,  increasing  slightly  with  isotopic  mass.  In  contrast  to  what  is 
expected  for  carbon  erosion,  the  yields  for  H  and  H2  ejection  are 
nearly  independent  of  the  mass  of  the  projectile.  Fig.  14  displays 
the  average  energies  of  the  ejected  particles  as  functions  of  the  im¬ 
pact  temperature  and  shows  that  they  are,  within  the  statistical 
uncertainty,  similar  for  all  three  impact  isotopes.  While  the  “re¬ 
flected”  energies  follow  the  impact  ones  almost  elastically,  the  en¬ 
ergy  of  the  sputtered  molecules  stays  constant,  as  in  case  of  beam 
irradiation  (Fig.  11). 


Fig.  13.  Isotope  specific  yields  of  H  (D,  T)  atoms  and  H2,  (D2,  T2)  molecules  upon 
irradiation  by  plasmas  composed  of  H  (green  squares),  D  (red  squares),  or  T  (blue 
squares)  at  various  temperatures.  Random  impacts  were  used  both  in  preparation 
of  the  surfaces  and  in  calculation  of  the  sputter  yield.  The  target  surface  was  kept  at 
a  temperature  of  300  K.  (For  interpretation  of  the  references  to  color  in  this  figure 
legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


Plasma  temperature  (K) 

Fig.  14.  Translational  energy  spectra  of  ejected  particles  upon  irradiation  as 
described  in  Fig.  13. 


3.  Outlook 

The  need  for  the  new  materials  in  fusion,  ITER  in  particular,  re¬ 
quires  studies  of  not  only  new  forms  of  carbon  (including  varia¬ 
tions  of  the  C  crystalline  and  polycrystalline  structures,  CEC,  and 
C  doped  with  boron  and  beryllium)  but  also  conceptually  new 
materials  like  tungsten  and  beryllium  and  their  compounds  and  al¬ 
loys  mutually  and  with  carbon  and  hydrogen.  Results  from  plasma 
irradiating  particles,  like  N(N2),  C,  W,  Be,  inert  gases,  as  well  as  iso¬ 
topic  effects  of  impacting  hydrogen  are  needed  for  a  more  com¬ 
plete  understanding  of  interactions  of  fusion  plasma  with 
plasma-facing  materials. 
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To  satisfy  these  needs,  development  of  new  potentials  that 
would  include  C,  W,  Be,  and  H  is  a  basic  prerequisite.  Further  devel¬ 
opment  of  the  hydrocarbon  potentials  to  improve  the  reaction 
dynamics  and  chemistry  of  complex  hydrocarbons  and  doped  car¬ 
bon  is  also  needed.  The  fast  expansion  of  computation  capabilities 
toward  petascale  machines  might  also  provide  the  capability  for 
more  realistic  simulations  in  the  development  of  new  methods, 
including  a  direct  quantum-mechanical  computation  of  the  re¬ 
quired  potentials,  introducing  electronic  excitations  and  charge 
transfer  during  the  impact  cascade  into  the  MD  simulations. 
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